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ABSTRACT 

We have performed high-resolution hydrodynamic simulations of the Lya forest 
in a variety of popular cold dark matter dominated cosmologies, including a low den- 
sity and a vacuum-dominated model. The fluctuation amplitude of these models is 
chosen to match the observed abundance of galaxy clusters at low redshift. We as- 
sume that the intergalactic medium is photoionized and photoheated by a uniform 
UV-background with the required amplitude to give the observed mean hydrogen 
absorption. We produce simulated spectra, analyze them by fitting Voigt profiles and 
compare line statistics with those obtained from high resolution observations. All mod- 
els give column density distributions in good agreement with observations. However, 
the distribution of line widths (the 6-parameter distribution) reflects differences in the 
temperature of the intergalactic medium between the models, with colder models pro- 
ducing more narrow lines. All the models with a low baryon density, f^ft 2 = 0.0125, 
are too cold to produce a 6-parameter distribution in agreement with observations. 
Models with a higher baryon density, fli,h 2 = 0.025, are hotter and provide better fits. 
Peculiar velocities contribute significantly to the line widths in models with low matter 
density, and this improves the agreement with observations further. We briefly discuss 
alternative mechanisms for reconciling the simulations with the observed 6-parameter 
distributions. 
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1 INTRODUCTION 

Neutral hydrogen in the intergalactic medium produces a 
forest of Lya absorption lines blueward of the Lya emis- 
sion line in quasar spectra (Bahcall & Salpeter 1965, Gunn 
& Peterson 1965). Hydrodynamic simulations of hierarchi- 
cal structure formation in a cold dark matter (CDM) dom- 
inated universe have been shown to be remarkably success- 
ful in reproducing a variety of properties of the observed 
forest (Cen et al. 1994, Zhang, Anninos & Norman 1995, 
Miralda-Escude et al. 1996, Hernquist et al. 1996, Wads- 
ley & Bond 1996, Zhang et al. 1997, Theuns et al. 1998ab, 
Dave et al. 1998, Bryan et al. 1998). The simulations show 
that the weaker Lya lines (neutral hydrogen column density 
iVjjj < 10 14 cm -2 ) are predominantly produced in the fila- 
mentary and sheet-like structures that form naturally in this 
structure formation scenario. Higher column density lines 
occur where the line of sight intersects a halo. 

These models assume that the intergalactic medium 
(IGM) is photoionized and photoheated by UV-light from 



quasars. The characteristic break in the rate of evolution of 
the number of lines below a redshift ~ 1.7, as observed by 
the Hubble Space Telescope, can then be explained by the 
decrease in the intensity of this ionizing background, itself a 
consequence of the rapid decline in quasar numbers towards 
lower redshifts (Theuns, Leonard & Efstathiou 1998, Dave 
et al. 1998). 

While the first simulations showed good agreement with 
the observed line statistics, more detailed studies at higher 
numerical resolution in a standard CDM universe produced 
a larger fraction of narrow lines than are observed (Theuns 
et al. 1998b, Bryan et al. 1998). Theuns et al. suggested that 
an increase in temperature of the IGM might broaden the 
absorption lines sufficiently to improve the agreement with 
observations. A higher gas temperature leads to more ther- 
mal broadening and in addition increases the Jeans mass. 
At low redshifts, the temperature of the gas responsible 
for the majority of Lya lines is determined by the balance 
between adiabatic cooling and photoheating. This causes 
the gas temperature to be a function of density, and this 
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Table 1. Models simulated 





Name 




n A 


C8 


n b h 2 


h 


r H , 


1 


S 


1 





0.50 


0.0125 


0.5 


HM/2 


2 





0.3 





0.85 


0.0125 


0.65 


HM/2 


3 


L 


0.3 


0.7 


0.90 


0.0125 


0.65 


HM/2 


4 


Sb 


1 





0.50 


0.025 


0.5 


HM 


5 


Ob 


0.3 





0.85 


0.025 


0.65 


HM 


6 


Lb 


0.3 


0.7 


0.90 


0.025 


0.65 


HM 



'equation of state' is well approximated by a power law, 
T = Tb(l + 5) 7_1 , where S is the gas overdensity (Hui & 
Gnedin 1997). To can be made higher by increasing the pho- 
toheating rate, or by heating the gas for longer by increas- 
ing the age of the universe. The first can be achieved by 
increasing the baryon density fi(,/i 2 , the latter by decreasing 
the matter density £l m . Using some simplifying assumptions, 
Hui and Gnedin (1997) obtain the scaling 

T c<[n i ,/i 2 /(n ro ft 2 ) 1/2 ] 1/1 - 7 . (1) 

Increasing Q,i,h 2 by a factor 2 and using fl m = 0.3 instead 
of 1 increases To by a factor > 2, which might be sufficient 
to obtain agreement with observations. 

In this Letter we use high resolution hydrodynamic sim- 
ulations to investigate the dependence of the 6-parameter 
distribution on cosmological parameters, for a given reion- 
ization history. This complements the work of Haehnelt & 
Steinmetz (1998), who quantified the dependence of the b- 
parameter distribution on the reionization history, for a crit- 
ical density Universe. 



2 SIMULATION 

We have simulated six different cosmological models, char- 
acterized by their total matter density fl m , the value of 
the cosmological constant Qa, the rms of mass fluctuations 
in spheres of 8/i -1 Mpc, erg, the baryon density Qi,h 2 and 
the present day value of the Hubble constant, Ho = 100ft 
km s~ Mpc - . The parameters for these models are sum- 
marized in Table |l|. Note that the Sb, Ob and Lb mod- 
els have a baryon density slightly higher than the value 
Q b h 2 = 0.0193 ± 0.0014, advocated by Buries & Tytler 
(1998) from measurements of the deuterium abundance in 
quasar spectra. Rauch et al. (1997) required a similarly high 
value Qbh 2 > 0.017 in order to make the measured flux 
decrement of the Lya -forest consistent with the observed 
intensity of ionizing photons from quasars. 

We model the evolution of a periodic, cubical region of 
the universe of comoving size 2.5/i _1 Mpc. The code used 
is based on a hierarchical P3M implementation (Couch- 
man 1991) for gravity and smoothed particle hydrodynamics 
(SPH, Lucy 1977, Gingold & Monaghan 1977, see e.g. Mon- 
aghan 1992 for a review) for hydrodynamics. It is a hy- 
brid between the HYDRA code of Couchman et al. (1995) 
and the APMSPH code, which were described in detail 
in Theuns et al. (1998b). In particular, we use the APM- 
SPH method for consistently computing SPH densities for 
particles in underdense regions. These simulations use 64 3 
particles of each species, so the SPH particle masses are 
1.65 x 10 6 (n 6 /i 2 /0.0125)(/i/0.5)" 3 M Q and the CDM parti- 
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Figure 1. Comparison between observed and simulated column 
density distributions at mean redshifts of 3 (top panel) and 2 (bot- 
tom panel). The data in the top panel are from Kim et al. (1997, 
hereafter KHCS, z = 3.35) and Hu et al. (1995, hereafter HKCSR, 
z = 2.9) and in the bottom panel from KHCS (z = 2.3) and from 
Petitjean et al. (1993, hereafter PWRCL, z = 2.0). The different 
models are indicated in the panel. 

cles are more massive by a factor ficDM/fit,- This resolu- 
tion is sufficient to simulate line widths reliably (Theuns 
et al. 1998b; note that numerical convergence will be even 
better in the hotter simulations). The simulations of Bryan 
et al. (1998) indicate that the absence of long waves will 
produce a small, but for our purposes unimportant, under- 
estimate of the widths of the absorption profiles. All our 
simulations were run with the same initial phases to mini- 
mize cosmic variance when comparing the different models. 

We assume that the IGM is ionized and photoheated 
by an imposed uniform background of UV-photons that 
originates from quasars, as computed by Haardt & Madau 
(1996). This flux is redshift dependent, due to the evolution 
of the quasar luminosity function. Haardt & Madau (1996) 
give two fits to the hydrogen ionization rate Phi as a function 
of redshift. We have used their qo = 0.5 fit for the critical 
density models and their qo = 0.1 fit for the open ones. This 
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Figure 2. Simulated ft-parameter distribution for models with low baryon density (top panels) and models with high baryon density 
(bottom panels) at redshifts 4, 3 and 2 (left to right). Superimposed with error bars are the observed data for z = 4 from Lu et at. (1996), 
HKCSR (z = 3) and KHCS (z = 2). We have only included lines with column density Nn t > 10 13 cm -2 that have an estimated error in 
b, berr < 0.256, for both simulations and data. 



amplitude of the flux is indicated as 'HM' in the Fhi column 
of Table |l|. In addition, for the low Qbh 2 models, we have 
divided the ionizing flux by two, indicated as 'HM/2' . The 
dependence on go reflects small differences in the complete- 
ness corrections of the observed quasar luminosity functions 
and is relatively unimportant for the simulations described 
here. The detailed expressions for the heating and cooling 
rates as a function of temperature and ionizing flux are taken 
from Cen (1992) with some minor modifications (Theuns et 
al. 1998b). Our analysis differs slightly from that in The- 
uns et al. in that we do not impose thermal equilibrium but 
solve the rate equations to track the abundances of Hi , Hll 
and Hel , Hell and Hem . We assume a helium abundance 
of Y = 0.24 by mass. 

We have used CMBFAST (Seljak & Zaldarriaga 1996) 
to compute the linear matter transfer function for each 
model. The amplitude of the transfer function is normal- 
ized to the observed abundance of galaxy clusters at z = 
using the fits erg = 0.52fi- - 46+o inm for fl A = and 



as = 0.52n-°- 52+013n ™ for fi A 
by Eke, Cole & Frenk (1996). 



1 — f2 m , as computed 



At several output times we compute simulated spectra 
along lines of sight through the simulation box. Each spec- 
trum is convolved with a Gaussian with full width at half 
maximum of FWHM = 8 km s _1 , then re-sampled onto pix- 
els of width 3 km s _1 to mimic the instrumental profile and 
characteristics of the HIRES spectrograph on the Keck tele- 
scope. We rescale the background flux in the analysis stage 
such that the mean effective optical depth at a given redshift 
in all models is the same as for the Ob model . The latter 
model has a mean absorption in good agreement with ob- 
servations (Ob has r c ff = 0.93, 0.33 and 0.14 at z=4, 3 and 
2) . Finally, we add to the flux in every pixel a Gaussian ran- 
dom signal with zero mean and standard deviation a = 0.02 
to mimic noise. We fit the continuum of the spectra using 
the method described in Theuns et al. (1998b). The absorp- 
tion features in these mock observations are then fitted with 
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Voigt profiles using an automated version of VPFIT (Car- 
swell et al. 1987). 



3 RESULTS AND DISCUSSION 

A comparison between the observed and simulated column 
density distribution functions (CDDFs) is shown in Fig. |l| 
The CDDFs for the different models are almost indistin- 
guishable and in addition are very similar to the observed 
distribution at both redshifts 2 and 3. In principle, the slope 
of the CDDF is sensitive to the amplitude of the mass fluc- 
tuations on small scales {e.g. Gnedin 1998). Figure |l| then 
shows that demanding consistency of the model with the ob- 
served abundance of galaxy clusters at z = automatically 
guarantees good agreement with the observed CDDF, for 
the CDM dominated cosmological models examined here. 

The distribution of line- widths for the different models 
is compared with observation in Fig. |^. We have selected 
those lines with column density 2Vhi > 10 13 cm -2 that have 
an estimated error from VPFIT b err < 0.256. For the data 
we have applied the same column density cut, and used the 
published uncertainties in 6 for each fitted Lya -line, where 
available. 

The models have different temperature-density rela- 
tions at z = 4, yet all of them provide fa-parameter dis- 
tributions which are in good agreement with the observed 
distribution at this redshift. However, at z — 3 there is a 
significant overproduction of narrow lines in all models with 
Qbh 2 = 0.0125, compared to the data. The models with the 
higher value Q,bh 2 = 0.025 agree much better with the ob- 
served fa-distribution, although the critical density Sb model 
still peaks at lower values than the observed distribution. 
The open models Ob and Lb, on the other hand, produce fa- 
parameter distributions which fit the observed distribution 
reasonably well. At redshift z = 2, the quality of the data 
is lower since the Keck spectrograph becomes less efficient 
in the blue. The high baryon models produce fa-parameter 
distributions which fit the observed distribution at z — 2 
better than the low ttbh 2 models. Again, the critical den- 
sity model peaks at significantly lower b than the observed 
distribution. 

We illustrate the main reasons for the cosmology de- 
pendence of the fa-parameter distribution in Fig. |^. We have 
recomputed simulated spectra for model S after imposing 
the equation of state of the significantly hotter model Ob. 
We did this by increasing the temperature of SPH parti- 
cles when they were cooler than the temperature given by 
the required equation of state. We then rescaled the flux so 
that the mean decrement was again the same as for the Ob 
model. Comparing this new model S-hot with S and Sb, we 
see that the main reason for the broader lines in Sb as com- 
pared to S is the higher temperature of the IGM in model 
Sb. However, the fa-distribution of model S-hot still peaks at 
lower values than that of model Ob. The reason for the re- 
maining differences is the contribution of peculiar velocities 
to the line-width for model Ob, as shown in Fig. Model 
'Ob-vel' is identical to model Ob, except that we have put 
all peculiar velocities to zero. This increases the fraction of 
narrow lines and the fa-distribution of this new model is al- 
most identical to that of models Sb and S-hot. The right 
panel in Fig. shows that neglecting peculiar velocities in 



model Sb decreases the fraction of narrow lines for this cos- 
mology, as shown by model Sb-vel. Consequently, the open 
model Ob has broader lines than the critical model Sb be- 
cause of the different way that peculiar velocities contribute 
to line broadening between the two models. The effects of 
peculiar velocity are similar in models Lb (not shown) and 
Ob. 

In summary, we have investigated the properties of the 
simulated Lya -forest in three popular cosmologies (criti- 
cal, open and vacuum-dominated) and compared simulated 
spectra analyzed in terms of Voigt profiles with observations. 
When the power spectrum is normalized to the abundance of 
galaxy clusters and the amplitude of the ionizing background 
is normalized to match the observed mean absorption, then 
all investigated models produce column density distributions 
which are well matched to the observations, at z — 2 and 
3. However, the distribution of fa-parameters is sensitive to 
the baryon density as well as to the matter density of the 
universe. Models with a higher baryon density Qbh 2 are hot- 
ter and hence the increased thermal broadening and Jeans 
smoothing shift the fa-parameter distribution towards larger 
b values, providing a better fit to the observations. Pecu- 
liar velocities contribute significantly to the broadening of 
absorption lines with column density Am ~ 10 13 cm -2 , in 
low matter density models. This broadening improves the 
agreement with observations further. In contrast, peculiar 
velocities narrow the lines in a critical density universe. The 
fa-distributions for our open and vacuum-dominated mod- 
els with Q m = 0.3 are almost identical. Of all our simu- 
lated models, the one that fits best has f^/i 2 = 0.025 and 
fl m — 0.3 and produces significantly more broad lines than 
the other models, yet there is a hint that the observed dis- 
tribution has an even larger fraction of broad lines. It is not 
clear whether the present data are sufficiently reliable that 
this is a real discrepancy. For example, some of the remain- 
ing differences may be due to differences in the Voigt profile 
fitting procedures between VPFIT and HKCSR or due to 
uncertainties in the continuum fitting (Rauch et al. 1993), 
although it seems unlikely that these effects would be large 
enough to reconcile our model with the data. If the discrep- 
ancy is real, then it is possible that feedback from e.g. star 
formation, not included in this investigation, might con- 
tribute to broadening the lines. Alternatively, a larger con- 
tribution from helium heating than assumed in our model 
might be required to boost the temperature further. 
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